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The filamentational instability of spatially broadband femtosecond optical pulses 
in air is investigated by means of a kinetic wave equation for spatially incoherent 
photons. An explicit expression for the spatial amplification rate is derived and 
analyzed. It is found that the spatial spectral broadening of the pulse can lead to 
stabilization of the filamentation instability. Thus, optical smoothing techniques 
could optimize current applications of ultra-short laser pulses, such as atmospheric 
remote sensing. 
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Recently, there has been a great deal of interes t 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 in investigating the nonlinear propa- 
gation of optical pulses in air. In order for the pulse propagation over a long distance, it is necessary 
to avoid filamentational instabilities that grow in space. Filamentation instabilities of optical pulses 
occur in nonlinear dispersive media, where the medium index of refraction depends on the pulse 
intensity. This happens in nonlinear optics (viz. a nonlinear Kerr medium) where a small modula- 
tion of the optical pulse amplitudes can grow in space due to the filamentation instability arising 
from the interplay between the medium nonlinearity and the pulse dispersion/diffraction. The fil- 
amentational instability is responsible for the break up of pulses into light pipes. It is, therefore, 
quite important to look for mechanisms that contribute to the nonlinear stability of optical pulses in 
nonlinear dispersive media. One possibility would be to use optical pulses that have finite spectral 
bandwidth, since the latter can significantly reduce the growth rate of the filamentation instabil- 
ity. Physically, this happens because of the distribution of the optical pulse intensity over a broad 
spectrum, which is unable to drive the filamentation instability with fuller efficiency, contrary to 
a coherent pulse which has a delta-function spectrum. In this Letter, we present for the first time 
a theoretical study of the filamentation instability of partially coherent optical pulses in air. We 
show that the spatial amplification rate of the filamentation instability is significantly reduced by 
means of spatial spectral broadening of optical pulses. The present results could be of significance 
in applications using ultra-short pulses for remote sensing of the near Earth atmosphere. 

The dynamics of coherent femtosecond optical pulses with a weak group velocity dispersion in 
air is governed by the modified nonlinear Schrodinger equatio n 4 ' 5 '?' 1 ^ 11 

a v/ + V 2 L v/ + /(|v/| 2 )v/ + /v|v/| 2 ^V = 0, (1) 

where i//(z,rjj is the spatial wave envelope, r_L = (x,y), and f(\ yr 2 \) = a\ yr\ 2 - e| v| 4 - Y\ ty\ 2K - 
Here a = 0.466, e = 7.3 x lO^cm 2 /^, y = 8.4 x lO^cm 2 ^- 1 )/^*^, and v = 1.2 x 
10~ 35 cm 2 ( K ~ 2 ^ /w^ K 2 ' for a pulse duration of 250 fs, and wq (in units of cm) is the beam waisfci^ 
(for a discussion of the approximations leading to Eq. (HJ), we refer to 9). We note that Eq. CQ) has 
been used in Ref. I 111 to analyze the multi-filamentation of optical beams. 

Following Ref. \V2[ we can derive a wave kinetic equation that governs the nonlinear propa- 
gation intense optical pulses which have a spectral broadening in space. Accordingly, we apply 
the Wigner-Moyal transform metho d 13 ' 14 ' 15 ' 16 . The multi-dimensional Wigner-Moyal transform, 
including the Klimontovich statistical average, is defined as 

p(z ' r± ' P) = (2^/ d2 Se ip -Hv*(z,r ± + a/2)Y(z,r ± -Z/2)), (2) 
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where p = (p x ,Py) represents the momenta of the quasiparticles and the angular bracket denotes 
the ensemble average 17 . The pulse intensity (| i//| 2 ) = I satisfies 



J d 2 pp(z,r^v). (3) 



Applying the transformation (2) on Eq. ©, we obtain the Wigner-Moyal kinetic equation 14 4 5 i 16 i 18 
for the evolution of the Wigner distribution function, 

d z p + 2p • V ± p + 2/(7) sin (± V± ■ V P ) P + 2vI K ~ l cos (j V± ■ V p ) p = 0. (4) 

Seeking the solution p = p(z,p) to Eq. ©, we may write p(z,p) = po(p)/(z), where po is 
an arbitrary function of p satisfying J d 2 pp = 1, and I(z) = I (2K - 2)l\2vl 2K ~ 2 z + (2K - 
2) 2K ~ 2 } 1 '( 1K ~ 2 > , with Iq = 1(0). Thus, the effect of a small but non-zero v is to introduce a 
slow fall-off in the intensity along the z-direction when K > 1 . Moreover, as V — > this solution 
reduces to / = To- 

We now consider spatial filamentation of a well defined optical pulses against small perturba- 
tions having the parallel wavenumber kn and the perpendicular wavevector kj_, by assuming that 
v is small so that k\\ ^> \d z \ for the background distribution. We let p = p(z,p) + p\(p) exp(iknz+ 
■ rj_) + c.c. and I = T(z) + 1\ exp(ihz + ik± ■ rj_) + c.c, where |pi | < p, \h \ <^ h and c.c. 
stands for the complex conjugate. We linearize (4) with respect to the perturbation variables and 
readily obtain the nonlinear dispersion equation 

, _ fr2 [f'(T)+iv(K-l)I K - 2 ]p(z,p-kj2) - [f>(I)-iv(K-l)F- 2 }p(z,p + kj2) 
J P k ll +2k ± -p-2ivI K - 1 ' K } 

which is valid for partially coherent femtosecond pulses in air. Here the prime denotes differenti- 
ation with respect to the background intensity I. 

We simplify the analysis by assuming that the perpendicular dependence in essence is one- 
dimensional. In the coherent case, i.e. p(z,p) = T(z)5(p — po), Eq. © yields 



*H = -2kp + iv(K+l)I K - l ±^k 2 [k 2 -2f(I)I\-v 2 (K-l) 2 P K - 2 , (6) 

where k represents the perpendicular wavenumber in the one-dimensional case. Letting h\ = 
—2kpo — iT in (6), where T is the filamentation instability growth rate, we thus obtair 19 



r = -v(K+ I)/*" 1 + Jk 2 [2f(I)I-k 2 ] + v 2 (K- l) 2 / 2 *- 2 , (7) 
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which reduces to the well known filamentation instability growth rate in a Kerr medium (i.e. v = 
and /(/) = al). We note that a nonzero v gives rise to an overall reduction of the growth rate. In 
Fig. 1 we have plotted a number of different curves for the growth rate in the coherent case. 

In the partially coherent case, we investigate the effects of spatial spectral broadening using the 
Lorentz distribution 

where A denotes the width of the distribution around the quasiparticle momenta pq. Inserting ([8]) 
into © and carrying out the integration in a straightforward manner, we obtain 



k\\ = -2kp + iv(K + l)I K - l +2ikA±^k 2 [k 2 -2f(T)I\ - v 2 (K-l) 2 P K - 2 . (9) 
With kn = —2kpo — tT the filamentation instability growth rate is 



r = - v (K + 1 ) I K - 1 - 2k A + J k 2 [2f (1)1 - k 2 } + v 2 (K - 1 ) 2/2K-2 . (io) 

In the limit A — > 0, Eq. © reduces to the dispersion relation ©, while for v = the dispersion 
relation © reduces to the standard expression for the filamentation instability growth rate 



r=-2kA + k^/2I f(I )-k 2 . (11) 

In Fig. 2 we have displayed the filamentation instability growth rate (TTOb . The effect of the finite 
width A of the quasiparticle distribution can clearly be seen. In particular, multi-photon absorption 
(here chosen to be a modest K = 3), determined by the coefficient v, as well as multi-photon 
ionization, represented by the coefficient y, combined with finite spectral width of the optical pulse 
give rise to a significant reduction of the filamentation instability growth rate. This is evident from 
Fig. 2, where the plotted normalized growth rate is reduced by as much as a factor of six, compared 
to the case of full coherence. 

In practice, optical smoothing techniques, such as the use of random phase plates 20 or other 
random phase techniques well suited for the results in the present Letter, have been used in inertial 
confinement fusion studies for quite some time (see, e.g. Ref. 1211) . Such spatial partial coherence 
controls are reproducible and can be tailored as to give a suitable broadband spectrum (as in, e.g. 



221 where optical vortices were generated). Thus, in the case of ultra-short pulse propagation in air, 
such random phase techniques can be used to experimentally prepare an ultra-short optical pulse 
for a long-distance propagation, and a large spatial bandwidth of optical pulses, in conjunction 
with multi-photon ionization and absorption, may drastically reduce (down to less than 20 % of 



the coherent value in the present study) the filamentation instability growth rate. This will lead 
to a greater long range stability, since the onset of strong optical pulse filamentation is delayed, 
resulting in several times longer stable propagation. A rough estimate based on the numbers found 
in the present Letter shows that an optical beam could propagate a distance as much as six times 
longer with proper random phasing. 

To summarize, we have investigated the filamentation instability of partially coherent fem- 
tosecond optical pulses in air. For this purpose, we introduced the Wigner-Moyal representation 
on the modified nonlinear Schrodinger equation and obtained a kinetic wave equation for optical 
pulses that have a spectral bandwidth in wavevector space. A perturbation analysis of the kinetic 
wave equation gives a nonlinear dispersion relation, which describes the filamentation instability 
(spatial amplification) of broadband optical pulses. Our results reveal that the latter would not be 
subjected to filamentation due to spectral pulse broadening. Hence, using partial spatial coherence 
effects for controlling the filamentational instability, femtosecond optical pulse propagation in air 
can be improved significantly. The result presented here is also indicative that optical smooth- 
ing techniques, as used in inertial confinement studies, could be very useful for ultra-short pulse 
propagation in air. This can help to optimize current applications of ultra-short laser pulses for 
atmospheric remote sensing over a long distance. 
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Fig. 1 . The coherent filamentation instability growth rate, given by Q, plotted for different param- 
eter values; all curves with Iq = 0.5, a = 1, and K = 3. The full thick line represents the standard 
filamentation instability growth rate for a nonlinear Schrodinger equation, i.e. v = £ = y = 0; the 
thin dashed curve has v = y = 0, while e = 0.5; the thin dotted curve has v = £ = and y = 0.5; 
the thin dashed-dotted curve has v = and £ = y = 0.5; the thick dashed curve has v = 0. 1 and 
e = y = 0; finally, the thick dashed-dotted curve has v = 0. 1 and £ = y = 1/2. 



Fig. 2. The partially coherent filamentation instability growth rate, given by (fTOl) . plotted for 
different parameter values; all curves with To = 0.75, a = 1, and K = 3. The full thick line again 
represents the standard filamentation instability growth rate for a nonlinear Schrodinger equation, 
i.e. A=v = £ = y=0; the thin full curve has v = £ = y = 0, while A = 0. 1 ; the thin dashed curve 
has £ = y = while v = 0.05 and A = 0. 1 ; the thin dotted curve has v = 0.05 and y = 0. 1 while 
£ = 0. The effects finite width of the background intensity distribution of the optical pulse, as well 
as the influence of the higher order nonlinearity and losses are clearly seen here. 
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